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£> Abstract 

Recent advances in attosecond spectroscopy techniques have fueled the interest in 
the theoretical description of electronic processes taking place in the subfemtosecond 
time scale. Here we study the coupled dynamic screening of a localized hole and a 
photoelectron emitted from a metal cluster using a semi-classical model. Electron 
density dynamics in the cluster is calculated with Time-Dependent Density Functional 
Theory and the motion of the photoemitted electron is described classically. We show 
that the dynamic screening of the hole by the cluster electrons affects the motion of 
the photoemitted electron. At the very beginning of its trajectory, the photoemitted 
electron interacts with the cluster electrons that pile up to screen the hole. Within our 
model, this gives rise to a significant reduction of the energy lost by the photoelectron. 



Thus, this is a velocity dependent effect that should be accounted for when calculating 
the average losses suffered by photoemitted electrons in metals. 
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1 Introduction 



Photoemission spectroscopy is one of the most important techniques used to study the 
structure of molecules, surfaces and solids [T]. It is based on the photoelectric effect, discov- 
ered more than 100 years ago by H. Hertz. Later, in 1905, Albert Einstein explained this 
effect as a quantum phenomenon, based on the emission of electrons from a target following 
the absorption of a photon ("quantum of light"). Photoemission spectroscopy has signifi- 
cantly contributed to the understanding of fundamental principles in solid state physics. 

In recent years, progress in laser technology has made possible the development of pho- 
toemission spectroscopy in the attosecond range (1 as= 10~ 18 s) [2J. Attosecond techniques 
permit access to the time scale of electron motion in atoms, molecules, and solids. Due to 
this experimental advance, there is a growing interest in the theoretical description of the 
dynamical electronic processes taking place in the subfemtosecond time scale [3j HJ 13 E]. 

In the present work, we study the electron dynamics during photoemission from small 
metallic clusters (or nanoparticles). Clusters represent a bridge between individual atoms 
and solid state materials [7]. Due to their large surface to bulk ratio, small metal clusters 
can exhibit rather unique features. For example, they frequently present interesting catalytic 
properties. Our choice of a finite size system as target in the photoemission process simplifies 
the theoretical analysis, but some of our conclusions are expected to remain valid in extended 
systems such as metal surfaces. 

We consider the case where one of the atoms in the metallic cluster undergoes core- 
electron photoemission. We focus our attention on the combined dynamic screening of a 
static localized core-hole and the photoemitted electron. We show that the presence of the 
hole left behind affects the many-body electronic dynamics in the cluster and therefore the 
emission dynamics of the photoelectron. For the description of the many-body response of the 
valence electrons in the cluster we use time- dependent density functional theory (TDDFT) - 
an ab-initio quantum-mechanical method. Our TDDFT methodology has already been suc- 
cessfully applied to study the dynamic screening of charges in finite-size systems [8j [9] and to 
the calculation of the energy transfer between particles and small gas-phase clusters [TOj E] • 
In the present calculations the motion of the photoemitted electron is described classically. 
This approximation is justified, provided that typical energies of the photons are in the 100 
eV range (XUV), resulting in relatively high energies of the photoemitted electrons, as the 
ones considered here. To analyze the role of the many-body screening effects we perform 
calculations using various approximations for the classical trajectory of the electron, includ- 
ing constant velocity studies and calculations with and without direct interaction between 
the ejected electron and the hole left behind. 

The problem we are addressing here has a long history in condensed matter physics. The 
dynamic relaxation of the Fermi sea after creation of a hole was analyzed in the context of 



X-ray photoemission by several authors [T2], [13] . Within the framework of linear response 
theory, Noguera et al. showed that the effective interaction between the core-hole and the 
photoemitted electron changes continuously from a statically screened potential for low- 
energy electrons to a completely unscreened potential for high-energy electrons [13]. They 
also showed that the double screening of hole and electron can occur with or without creation 
of plasmons according to the kinetic energy of the emitted electron. Here we go beyond linear 
theory in the description of the dynamic screening of charges in the photoemission process 
by using propagation of electronic wave-packets with TDDFT to compute the response of 
the valence electrons. 

2 Theory 



In the present study metallic clusters are described using a spherical jellium model (JM). 
Despite its simplicity, the JM can be very useful in the interpretation of photoemission data 
from metal clusters, as recently shown in Ref. [H]. In the JM, the core ions are substituted 
by an homogeneous background of positive charge with a density defined by 

n+(r)=n (r s )G(R d -r), (1) 

where R c i is the radius of the cluster, 9(x) is the Heaviside step-function and no(r s ) is the 
constant bulk density, which depends only on the Wigner-Seitz radius r s (1/uq = Airrl/3) 
|15j . The latter is the only parameter in the JM. The number of electrons in a neutral cluster 
is N = (R c i/r s ) . For simplicity, we only consider closed-shell clusters in our calculations. 

In order to obtain the ground state electronic density of the cluster n(r) we use the spin- 
restricted density functional theory (DFT) [16] and solve Kohn-Sham (KS) equations [TT] : 

l-^V 2 + V eff ([n],r)\ Vl = e m , (2) 

where Si are the eigenvalues of the KS equations and ifi are the one-electron wave functions. 
Please, notice that, unless otherwise specified, we use Hartree atomic units (a.u.) throughout 
the paper. The effective potential is composed of three terms: 

V eff ([n], r) = V ext (r) + V H (r) + V xc (r), (3) 

where V ext (r) is the external potential created by the positive background, Vuij) is the 
Hartree (or Coulomb) potential created by the electronic density, and V xc (r) is the exchange- 
correlation potential, calculated in our case in the local-density approximation (LDA) with 
the Perdew-Zunger parametrization of Ceperley- Alder exchange and correlation potential [18] . 
The electronic density n(r) is given by a sum over occupied wave functions: 



idocc 



2^ |^(r)| 2 , (4) 



where the factor 2 stands for spin degeneracy. Results for the ground state of small metal 
clusters have been thoroughly discussed in the literature [191 EDJ [21]. The effective KS 
potential and the ground state electronic density for one of the clusters considered here are 
shown in Figure [l|a). 



For the description of the photoemission process we use a semi-classical model. We 
consider fast photoemitted electrons, moving with velocities much higher than the Fermi 
velocity of the cluster electrons. Therefore, the movement of the photoelectron can be 
represented classically. In parallel, valence electron dynamics in the cluster is investigated 
by means of TDDFT (22] ■ In TDDFT, the evolution of the electronic density n(r) in response 
to the field of the moving electron is calculated solving the time-dependent KS equations: 



■ 9 , N 



-V 2 + V ext {v) + V H (r,t) + V xc (r,t) + AV(r,t) } ^(r,i) 



(5) 



where AU(r, t) is the change of the external potential due to the photoemission process (see 
below). The exchange-correlation potential V xc is calculated with the standard adiabatic 
local density approximation (ALDA) with the parametrization of Ref . [TB] . 

The time-evolving electronic density of the excited cluster is obtained from the time- 
dependent KS orbitals (fi(r,t), in a way similar to Eq. (HI). The time-dependent KS wave 
functions are obtained by propagating the initial wave functions (fi(r,t ) = (pi(r) using the 
split-operator technique. Due to the presence of the photoemitted electron the problem loses 
its spherical symmetry and the use of cylindrical coordinates (p, z) becomes necessary. A 
detailed description of the numerical procedure can be found in Refs. [231 EH EE] • 

In our model, we do not consider explicitly the interaction with an external electromag- 
netic field. Thus, an electron with high kinetic energy and a static hole are created at t — 
at the center of the spherical cluster. The scheme of the process is shown in Figure [lib). 
The photoemitted electron (el) is modeled as a negative point charge that moves along the 
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Figure 1: a) Upper panel, ground state electronic density in units of the positive back- 
ground density (dashed line, r s =4) for a cluster containing 20 electrons. Lower panel, the 
corresponding effective Kohn-Sham potential and occupied energy levels, b) Sketch of the 
photoemission process. An electron and a hole are created at the center of the spherical 
cluster at t = 0. Both are represented by classical point particle and the electron starts to 
move along the z- axis with velocity Vq. 



2-axis [p = 0, z e i(t)]. It is worth noting that this photoemitted electron is not one of the 
cluster valence electrons, but an extra electron coming from the inner shell of a hypothetical 



atom sited in the center of the cluster. The electron motion is calculated in two different 
approximations as described in detail below. The hole (h) is represented by a positive point 
charge at a fixed position [p = 0, Zh — 0]. The potential created by these charges and acting 
on the rest of the electrons in the cluster is AV = V e i + Vh, where 

Vd = Tt — 7T\ \2 i 2ii/2 (5) ( t )' ( 6 ) 

[{zei(t) - z) 1 + p l \ v l l 

and 

To address the effect of the many-body dynamics in the cluster on the energy loss expe- 
rienced by the ejected electron in well-defined conditions, we first study a simplified case in 
which the photoemitted electron is assumed to move with constant velocity v, i.e., z e i(t) = vt. 
This allows us to isolate the effects related to the dynamics of the screening processes from 
other possible effects associated with the details of the trajectory. Here, v corresponds to the 
final velocity of the electron if the photoemission process would take place in vacuum, which 
is considered as a good approximation for the average electron velocity during its movement 
through the cluster. Thus, the direct interaction between the electron and the hole is not 
taken into account in this case. However, as we will see below, the screening of the hole still 
has an important influence on the energy loss by the photoemitted electron. In this case the 
energy loss is calculated from the integral 



POO 

E loas = -v / F?'{t) dt, 

Jo 



where 



F z (t) = 2vr J dpdz P y— -^— — -gp [z el (t) - z] (9) 

is the ^-component of the force created by the cluster on the emitted electron. It is impor- 
tant to note that Ei oss includes the energy necessary to eject the electron from the cluster 
(an adiabatic contribution), as well as non-adiabatic contributions due to the creation of 
electronic excitations in the cluster during the emission process. 

In a second step, we perform a more refined treatment in which the direct electron-hole 
interaction is included and the trajectory z e i(t) is calculated using the classical equations of 
motion: 

dzd/dt = v(t), z el (t = 0) = (10) 

dv/dt = F 2 io4 (t), v{t = 0) = v Q . (11) 

In Eqs.(10) and (11), F l ° l is the total force felt by the moving electron: 



?tot/j.\ T^clsu\ z el\t) 



™ = ™ - M t)" + ^ - (12) 

The first term corresponds to the interaction with the cluster given by Eq.pl). The second 
term stands for the force due to the interaction between photoemitted electron and the 
core hole left behind. The electron-hole interaction in our study is given by the regularized 
Coulomb potential V e i-h- 

Vd ' k = " I l 2 - (13) 

We use a 2 = 0.5 to avoid divergence at time t — 0. 



3 Results and discussion 



In what follows, we center our discussion on the force experienced by the photoemitted 
electron due to the interaction with the cluster. This force is given by Eq.pl) and allows 
to study important aspects of the electron density dynamics in the cluster. We study such 
force in the two approximations mentioned in the previous section for the motion of the 
(classical) photoemitted electron. In the first approximation, in which the electron moves 
with constant velocity, we can conveniently identify the effect of the hole screening on the 
movement of the photoemitted electron. To quantify the effect of the dynamic screening 
of the hole we calculate the work performed by the force Ff s along the electron trajectory. 
This quantity is directly linked with the energy loss of the ejected particle. In the second 
approximation, the velocity of the electron is allowed to vary according to Newton laws. 
This approximation might be closer to the real photoemission process. Also in this case we 
find an important influence of the hole-screening dynamics on the force experienced by the 
emitted electron and, thus, on the energy loss during the photoemission process. 



3.1 Constant velocity approximation 

In the constant velocity approximation, we calculate the cluster induced force in two 
different cases, namely with a localized hole at the center of the cluster (potentials in Eqs. 
(6) and (J7|) are included in the calculations) and without the hole (only the potential in Eq. 
(6) is included). In this approximation the direct interaction between hole and electron is 
not included. In spite of this, we find that the presence of the hole modifies the electron 
dynamics in the cluster because there are two different charges to be screened. We have 
studied four Na clusters (r s = 4) with 20, 58, 106 and 556 electrons. We consider three 
different velocities of the photoemitted electron: 1 a.u, 1.5 a.u. and 2.5 a.u. Figure [2] shows 
the cluster induced force acting on the photoemitted electron as a function of the electron 
position in the cases v = 1 a.u. and v = 1.5 a.u. for the small cluster with 20 electrons. 
From the present results it follows that, in the case of the presence of the hole, the cluster 
induced force on the photoemitted electron has positive value at short times for the two 
chosen velocities. This indicates that the cluster response tends to accelerate the electron 
at the very beginning of its movement. The repulsive cluster induced force is related to the 
screening of the hole. More precisely, at the beginning of its movement, when close to the 
hole, the electron is repelled from the hole vicinity by the cluster electrons that arrive to 
screen the hole. When the hole is not included in the calculation, the above effect is not 
observed and the emitted electron is mainly decelerated all along its trajectory [26]. This 
deceleration is due to two effects. First, within the cluster, the electron suffers the stopping 
characteristic of any charged particle moving in an electron gas [11| . Second, as the electron 
approaches the cluster surface, we can clearly see the contribution to the forces associated 
with overcoming the surface potential step. 

These two decelerating contributions are difficult to disentangle for very small clusters, 
like those in the main panels of Figure [2j However, the force experienced by a particle moving 
inside a large jellium cluster reaches a stationary regime and oscillates around a mean value. 
This can be seen in the inset of Figure pq (a) for a cluster containing 556 electrons. The mean 
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Figure 2: Cluster induced force [Eq. (J9l)] acting on the electron moving away from the center 
of the cluster (N=20, r s =4) as a function of the electron position. We consider the cases 
without and with a localized hole at the center of the cluster. The electron moves with a 
constant velocity: a) v = 1 a.u. (Inset: Cluster induced force for a larger cluster with 556 
electrons, the hole is not included. The horizontal dotted line corresponds to the stopping 
power in an homogeneous electron gas with r s =4.) and b) v — 1.5 a.u. Vertical dashed lines 
correspond to the radius of the cluster, R d = 10.86 a.u. 



value of the force is the so-called stopping power and only depends on the electron density. 
In our case, for r s =4, is around 0.055 a.u. for a negatively charged particle moving with a 
velocity of 1 a.u. [TO] . 

The influence of the hole screening on the moving electron can be better analyzed if we 

Here, 



look at the difference of forces shown in the Figure pi F^ 



Fteli 



pels I 



rpcls 
r h,el 



(z) (F°j s ) is the cluster induced force on the photoemitted electron calculated with 



(without) explicit inclusion of the positive point charge at the center of the cluster. With 
this definition, F£ ls (z) is the force felt by the photoemitted electron due specifically to the 
cloud of electronic density that dynamically screens the hole. This quantity is shown in 
Figure [3] for all clusters considered and for two different velocities. 

As it is seen from the graphs, the effect of the screening of the hole is larger in the case 
of the smaller electron velocity 1 a.u. - Figure [3ta). This is related to the time that the 
photoelectron spends in the neighborhood of the hole and to the characteristic time of the 
hole screening. The slow photoelectron stays near the hole long enough for the screening 
of the hole to be performed. Therefore, it experiences a large force due to the piling up 
of electronic charge around the hole. The fast electron, however, leaves the hole at short 
times, which are not enough for a significant piling up of screening charge. Hence, the 
effect associated with the hole screening becomes smaller the higher the electron velocity. 



It is worth noting that, for the slow electron, Fj 



els I 



is almost identical for the two largest 



clusters considered here and it is very small for z e i > 25 a.u. Both observations indicate 
that the screening of the hole is well established and basically reaches its stationary value 
at the corresponding time scale. For the faster electron, however, the value of F£ ls (z) at 
large z e \ is different for different cluster radii. This is linked to the time-evolution of the 
screening density which still goes on by the time the electron reaches the cluster boundary. 
These conclusions are corroborated by the induced electron density dynamics plots (Figure [5]) 
discussed below. 
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Figure 3: Component of the cluster induced force acting on the moving electron due to 



the dynamic screening of the hole by the electronic charge of the cluster, Fj 



els { 
h \ 



. Results 

are shown as a function of the electron position for jellium clusters (r s =4) of different size 
comprising N = 20, 58, 106, and 556 electrons. Electron velocities are: a) v — 1 a.u, b) 
v = 2.5 a.u. Vertical dashed lines correspond to the radii of the clusters. 



In order to quantify the effect of the hole screening, we calculate the electron energy 
loss due to the interaction with the cluster electrons for the case of N = 20 and for the 
velocities 1 a.u, 1.5 a.u., and 2.5 a.u. of the photoemitted electron. The energy loss Ei oss , 
given by Eq. (J8|, is defined as the work performed by the cluster induced force acting on 
the moving electron. The results are summarized in Table [TJ in which the difference in the 
cluster induced energy loss AEi oss , with and without the hole, is also shown. 



Table 1: Estimated energy loss [Eq. rt8j)] for a photoemitted electron as a function of its 
velocity. The electron is assumed to follow a trajectory with constant velocity. The difference 
in energy loss AEi oss for the case of absence and presence of a hole in the center of the cluster 
is also given. 



N=20 



Ei oss , a.u. 



without hole 
with hole 



AE lc 



a.u. 



v — 1 a.u. 



0.60 
T3T 



0.23 



v 



1.5 a.u. 



2.5 a.u. 



0.44 
1)W 



0.22 
"016" 



0.14 



0.06 



The presence of the hole reduces the cluster induced energy loss for all velocities. The 
value of AEi oss also shows that the effect of the hole screening is more significant the slower 
the electron. The energy loss of the electron moving at 1 a.u. decreases almost by a factor of 2 
when we include the hole screening in the process. An interesting consequence is that, at low 
velocities, the effects associated with the hole screening might become crucial to determine 
if the photoemission process can indeed take place or not. For example, the kinetic energy 
of the slowest electron considered in Table [I] is 0.5 a.u.. Since the energy loss in the case 
without hole is 0.6 a.u., this electron cannot be photoemitted from the cluster. However, in 
the presence of the hole the photoemission becomes possible. 

The study of the electron dynamics during photoemission in the constant velocity ap- 
proximation leads us to two conclusions: 1) the screening of the hole by the cluster electrons 



8 



leads to a repulsive (accelerating) force acting on the photoemitted electron at the beginning 
of its movement; 2) the effect of the hole screening is reduced for faster (more energetic) 
photoemitted electrons. 



3.2 Varying velocity approximation 



The results discussed so far are obtained using a simple model in which the photoemitted 
electron moves with a constant velocity. In a real photoemission process, however, the 
velocity varies due to the different elastic and inelastic forces acting on the electron. In order 
to be sure that none of the effects discussed above is an artifact of the model and to prove our 
conclusions, we simulate the photoemission process in a more realistic second approximation. 
In this approximation, the velocity and coordinate of the electron are dependent on time, 
according to equations (10), (11), and (12). Electron and hole interact via a regularized 



Coulomb potential (13). 
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Figure 4: Cluster induced force [Eq. d9J)] acting on the electron moving away from the center 
of the cluster (N=20, r s = 4), as a function of the electronic position in three cases: 1) There 
is no hole at the center of the cluster and the initial velocity of the electron is 1.5 a.u. (green 
line); 2) There is a hole at the center of the cluster, but the photoemitted electron does not 
interact directly with the hole. The initial velocity of the electron is 1.5 a.u. (red line); 3) 



There is a hole at the center of the cluster and a direct electron-hole interaction (Eq. 13 ) is 



included. The initial velocity of the electron is 2.25 a.u. (red line). The vertical dashed line 
shows the radius of the cluster, R c i = 10.86 a.u. 



In Figure [4] we show the results for the small cluster with 20 electrons in this varying 
velocity approximation. The force felt by the moving electron due to the interaction with the 
cluster electrons [Eq. ^} is calculated for three different cases. In the first case, we do not 
include the hole in the cluster. In the second case, there is a hole to be screened, but we do 
not include the direct interaction between the hole and the photoelectron when performing 
the trajectory calculation. In the third case, the direct interaction between photoelectron 
and hole is included. In the first two situations, the initial velocity of the electron is set to 



1.5 a.u, while in the third situation is set to 2.25 a.u. This difference in velocity corresponds 
to the energy of the electron-hole interaction. Indeed, with the parameter o? = 0.5 a.u. 



used here in Eq.(13), the "binding energy" [27] of the electron is around 1.4 a.u. Taking 
into account this binding energy, the initial velocity 2.25 a.u. for an electron photoemitted 
in vacuum and interacting with the hole leads to a final velocity ~ 1.5 a.u. Therefore, our 
choice of initial velocities allows a direct comparison of the two cases, with and without 
direct electron-hole interaction in the case of photoemission inside the cluster. 

One can see from Figure [4] that the behavior of the cluster induced force in this more 
realistic model is similar to that of the simple model considered before (Figure [2]). Whenever 
the hole screening is taken into account, there is an acceleration force acting on the electron 
at the beginning of its trip. Similarly to the constant velocity approximation, performing cal- 
culations along a more realistic trajectory with different launch velocities, we also found that 
the effect of the hole screening decreases when the initial velocity of the electron increases. 

Moreover, we found that the cluster induced force is very similar independently on 
whether electron and hole directly interact with each other or not. However, this is valid 
only if the final velocity of the photoemitted electron interacting with the hole is equal to 
the final velocity of the electron not interacting with the hole. This shows that, as far as the 
final energy of the photoemitted electron is the same, the cluster induced force acting on the 
photoemitted electron is mainly affected by the presence or absence of the hole screening 
and not by the details of the electron trajectory nearby the hole. 

3.3 Time evolution of electronic density 

Continuing the discussion on the electron density dynamics in the cluster, we illustrate 
the effect of the coupling of both processes - dynamic screening of the hole and dynamic 
screening of the moving electron. Figure [5] shows the time evolution of the electronic density 
of the spherical cluster with N = 106 electrons. The hole and the electron are created at 
time t — at the center of the cluster (-2=0) and the electron is moving along the positive 
part of the z-axis with a constant velocity v = 1 a.u. The induced electronic density close 
to the z-axis An = n(po,z,t) — n(po,z, 0), where po = 0.02 a.u, is plotted in units of the 
background density n . 

Panel (a) shows the results of a direct TDDFT calculation of the induced electronic 
density for the cluster with a static hole at the center and a photoelectron moving away 
from the center of the cluster along the z-axis. Panel (b) shows the induced electronic 
density obtained as a sum of two different contributions: 

An{t) = 6n h {t) + 8n d {t). (14) 

Here, 5rih(t) is the TDDFT result for the induced electron density due to the appearance of 
only a localized hole at the center of the cluster. Similarly, Sn e i(t) is the TDDFT result for 
the induced electron density in response to a photoelectron moving from its center, where no 
hole is present. Therefore, panel (b) shows a linear superposition of the electronic charges 
screening the static hole and the moving photoelectron. In the inset of both graphs we show 
the time evolution of the electronic density at a given point (p=0.02 a.u., z=0.2 a.u.). 

The white area in the main plots shows a depletion of the electronic density in the cluster 
that roughly follows the trajectory of the electron. It is due to the Coulomb repulsion between 
the moving electron and the rest of the electrons in the cluster. Black arrows indicate the 
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Figure 5: Time evolution of the electronic density of the spherical cluster in which the hole 
and the electron are created at time t = at the center of the cluster and the electron moves 
away along the positive part of the z-axis. The induced electronic density is shown close to the 
symmetry z-axis (p=0.02 a.u.,z). The time evolves along the vertical axis. Color map shows 
the change in density in units of the background density n . The color scale is limited to a 
maximum value of 20 in order to reveal the effects in the regions where the induced density 
is small. The induced density above this value is shown in green. The actual maximum 
value of the induced density is about 46 — 50 in units of background density. It corresponds 
to the small z region around the position of the hole. Panel (a) shows the results of the 
TDDFT calculation of the complete system. In panel (b) the induced density is calculated 
as a sum of two contributions (see the text for the explanation). Cluster parameters are 
r s = 4, N = 106, R c i = 18.93 a.u. The velocity of the electron is constant and is equal to 
v = 1 a.u. Insets: Profile of the plot along the time axis at (p=0.02 a.u., z=0.2 a.u.). 

time at which the screening of the hole is fully developed, i.e., the induced electron density 
in the close vicinity of the hole roughly integrates to one. This time is also shown in the inset 
of each plot and is equal to 11 a.u. and 8 a.u. for the cases (a) and (b), respectively. Thus, 
there is a delay in the TDDFT screening of the hole as compared to the linear superposition 
case. Moreover, comparing the charge distribution for negative and positive values of z, we 
can see a clear asymmetry in the screening charge for the TDDFT calculation with both hole 
and electron simultaneously included. This asymmetry is absent in panel (b), corresponding 
to the linear superposition of electron and hole separate screenings, and clearly indicates 
that the dynamics of the hole screening is affected by the presence and movement of the 
emitted electron. Therefore, we can conclude that the TDDFT calculation, considering both 
the hole at the center and the electron photoemitted from the center of the cluster, includes 
a combined effect of the dynamic screening of both particles in the relaxation processes in 
the cluster. This combined effect is also visible in the oscillations of the electronic density, 
where the periods of these oscillations are slightly different for the two cases considered. 

4 Conclusions 



In this study a semi-classical model was used to describe the dynamic screening of a mov- 
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ing photoelectron and that of the localized core-hole left behind as a result of interaction of 
XUV pulses with small metal clusters. The motion of the photoemitted electron is described 
classically and the electron dynamics in the clusters is studied using the Time-Dependent 
Density Functional Theory, TDDFT. 

We have shown that, when the hole is explicitly included in the calculation, the pho- 
toemitted electron is accelerated by the cluster electrons that pile up nearby the cluster 
center to dynamically screen the hole. This effect is observed by comparing the forces act- 
ing on the photoemitted electron due to the interaction with the cluster in which a hole 
is present or absent at the center. In order to quantify the effect of the hole screening we 
have calculated the energy loss of the photoelectron. We have shown that the presence of 
the hole reduces significantly the cluster induced energy loss and that this effect is velocity 
dependent. The higher is the energy of the photoemitted electron the smaller is the effect 
induced by the hole screening. 

These conclusions were obtained using a relatively simple approximation in which the 
photoemitted electron moves with constant velocity. The conclusions are proven to remain 
valid when the interaction between photoemitted electron and core-hole left behind is in- 
cluded in the calculation and the velocity of the electron is allowed to vary with time. We 
have illustrated the time evolution of the electron density in the cluster during the photoe- 
mission process and we have shown that the TDDFT calculation allows us to see the coupled 
effect of the screening of both the hole and the electron in the relaxation processes inside 
the cluster. 

The semi-classical model used here allows for a detailed analysis of the effect of the dy- 
namic screening of the hole. However, the simplicity of the model and the classical treatment 
of the photoemitted electron also prevents a direct translation of our findings to the experi- 
mental situation. Thus, a clear understanding of the implications of the present results for 
photoemission experiments is still an open question that requires further work. 
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